Direct mapping of the formation of a persistent spin helix 
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The spin-orbit interaction (SOI) in zincblende semiconductor quantum wells can be set to a 
symmetry point, in which spin decay is strongly suppressed for a helical spin mode. Signatures of 
such a persistent spin helix (PSH) have been probed using the transient spin grating technique, 
but it has not yet been possible to observe the formation and the helical nature of a PSH. Here we 
directly map the diffusive evolution of a local spin excitation into a helical spin mode by a time- 
and spatially resolved magneto-optical Kerr rotation technique. Depending on its in-plane direction, 
an external magnetic field interacts differently with the spin mode and either highlights its helical 
nature or destroys the SU(2) symmetry of the SOI and thus decreases the spin lifetime. All relevant 
SOI parameters are experimentally determined and confirmed with a numerical simulation of spin 
diffusion in the presence of SOI. 



Conduction-band electrons in semiconductors experi- 
ence SOI from intrinsic pQ and extrinsic sources, lead- 
ing to spin dephasing, current-induced spin polarization 
and spin Hall effects [2]. These physical mechanisms are 
of great fundamental and technological interest, recently 
also in the context of topolocial insulators [3] and Majo- 
rana fermions [IJ|5]. Intrinsic SOI arises from an inver- 
sion asymmetry of the bulk crystal (Dresselhaus term) 
and of the grown layer structure (Rashba term). In a 
quantum well (QW), these two components can be tai- 
lored by means of the confinement potential [6 , and the 
Rashba SOI can be externally tuned by using gate elec- 
trodes [TIE]. In general, SOI leads to precession of elec- 
tron spins. In the diffusive limit, in which the scattering 
length is much smaller than the spin-orbit (SO) length 
Ago? a random walk of the spins on the Bloch sphere will 
dephase a non-equilibrium spin polarization [9]. 

Of special interest is the situation in a two-dimensional 
electron gas (2DEG) with balanced Rashba and Dressel- 
haus contributions [T0HT3] . There, the SOI attains 
SU(2) symmetry and the spin polarization of a helical 
mode is preserved. The reason for this conservation of 
the spin polarization is a unidirectional effective SO mag- 
netic field Bso, which depends linearly on the component 
of the electron momentum along a specific in-plane di- 
rection. This causes the precession angle of a moving 
electron to vary linearly with the distance traveled along 
that direction, irrespective of whether the electron path 
is ballistic or diffusive [lOj [11]. In such a situation, a 
local spin excitation is predicted to evolve into a helical 
spin mode termed PSH [Fig. [TJa)]. Transient spin grat- 
ing measurements [6] showed that a spin excitation with 
a spatially modulated out-of-plane spin component de- 
cays with two characteristic lifetimes that correspond to 
two superposed spin modes of opposite helicity. 



Here we directly measure the diffusive evolution of a 
local spin excitation into a PSH by time-resolved Kerr 
rotation microscopy [Fig. [ljb)]. We employ a pump- 
probe approach, in which a circularly polarized pump 
pulse excites electrons into the conduction band of a 
(OOl)-grown GaAs/AlGaAs QW with their spins polar- 
ized along ^ 1 1 [001] . The out-of-plane spin polarization S z 
is then measured by a probe pulse delayed by a time t, 
using the polar magneto-optical Kerr effect. The posi- 
tion of the incident pump beam is scanned to record the 
spatial spin distribution S z (x,y) at time t. We define x 
along the [110] and y along the [110] direction. The SOI 
of the 2DEG is tuned close to the SU(2) symmetry point, 
M ~ — by controlling the Rashba (a), the linear 
Dresselhaus and the cubic Dresselhaus (^3) SO cou- 
pling coefficient via asymmetric modulation doping on 
the two sides of the QW. 

The experimental observation of the PSH is exempli- 
fied by three maps of S z (x,y) recorded at different t 
[Fig. [IJc)]. The first map at t = 10 ps still shows the lo- 
cal excitation of S z > centered at x = y = 0. Because 
of the initially rapid spin diffusion, the Gaussian shape 
of S z (x,y) is already broader than the size of the fo- 
cused pump-laser spot. Spins further diffuse in the (x, y) 
plane, but the second and the third map (recorded at 
t = 240 and 840 ps) in addition feature alternating stripes 
of S z (x,y) > and S z (x,y) < caused by spin preces- 
sion about Bso- To explain this unidirectional oscillation 
along the ^-direction, Bso must be more strongly corre- 
lated with k y than with k x (k x and k y are the components 
of the electron wave vector k). With our definition of a 
and j3 and from the symmetry of Bso (see supplemen- 
tary information), it follows that a and (/3i — /3 s) must 
have the same signs and that therefore the x-component 
of Bso is much larger than the ^/-component. For op- 
posite signs of a and (j3\ — ^3), the PSH would oscillate 
along the x-direction and the ^/-component of Bso would 
be larger. 
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The formation of the PSH is best illustrated if spin 
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FIG. 1: Direct mapping of the persistent spin helix formation, (a), Diffusive expansion of a local spin excitation 
(top), where the spin polarization evolves into a PSH mode (bottom). The arrows and the colors indicate the direction of S 
and the magnitude of S z respectively, (b), Schematic of time-resolved Kerr rotation microscopy (see methods) and profile of 
conduction-band (CB) energy of the 12-nm-wide GaAs/AlGaAs QW sample investigated, (c), Experimental observation of 
the PSH. Spatial maps of S z (x,y) are shown for three different times t. 



dynamics is tracked in space and time. For that pur- 
pose, we position the pump pulse at x = and scan 
it along the ^-direction. A collection of such line scans 
S z (y,t) recorded at various t is shown as color-scale plot 
in Fig. [2|c). Starting with S || z, the excited spin dis- 
tribution expands along ±y, and thereby S z starts to 
oscillate with y. As we will discuss in the following, this 
oscillation is indeed the footprint of a helical spin mode. 

A cartoon of helical spin modes is shown in Fig. [2ja) . 
The spin polarization rotates by an angle that depends 
linearly on the position with the direction of rotation 
defining a helicity (cj + or uj~). Therefore, not only S z 
but also S y oscillates with y [Fig. [2|a)]. The helicity of 
the emerging spin mode depends on the sign of the cross- 
correlation (Bso yX k y } : i.e., on whether Bso, x is positive 
or negative for k y > 0, and is given by the absolute sign of 
a + Pi — Ps (see supplementary information). The helic- 
ity cannot be directly determined from S z (y), as for both 



modes, S z oscillates in the same way [Fig. [2f a)]. This is 
especially true for experiments employing the transient 
spin-grating method, in which as an additional compli- 
cation, spin waves with both helicities are excited simul- 
taneously [6]. 

We uncover the helical nature of the measured spin 
mode by rotating in-plane spin components S y out of 
plane with the help of an external magnetic field B || x 
[Fig, fflb)]. The effect of B x on S z (y,t) is shown in 
Fig. [2^d). S z (y,t) still oscillates with y - indicating that 
the PSH is preserved - but the phase of the oscillation 
now shifts with t (dashed line). This shift can be un- 
derstood from the spin precession about B X1 as shown 
in Fig. [2^b): For an cj + mode and B x < 0, the posi- 
tion yo of equal phase shifts toward — y with increasing 
t. As we observe the same direction in the measurement 
[Fig. J2^d)] , we conclude that our sample supports an uo + 
mode, which means that a + j3\ — Ps must be positive. 
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FIG. 2: Helical spin modes and the PSH. (a), The spins of the two helical spin modes cj + and u~ rotate with opposite 
helicity as the electron moves along y. The S z (y) of the two modes is the same, (b), Spins precess about an external magnetic 
field B x such that the cj + mode is shifted towards — y if B x < is applied, (c), (d), Maps of S z (y,t) show the time evolution 
of the PSH for B x = (c) and B x = — IT (d). S z (y) oscillates in both cases, but for B x — —IT, the position of constant 
phase, 2/0, is shifted with t (marked by the dashed line). 



With this, we have directly determined the sign of the 
cross-correlation (Bso yX k y ). 

Figure [2jd) also shows that for all positions y, S z (y,t) 
oscillates in t with the same frequency v. Therefore, 
S z (y,t) represents a collective precession of the uj + mode 
about B x . According to v = \gjiiBB x /h\^ we determine 
the electron g- factor g to equal —0.17 (based on the QW 
thickness, we assume g < 0; hb is the Bohr magneton, 
and h = 2tyH Planck's constant). 

We now quantitatively analyze S z (y,i) to extract in- 
formation on the SOI and on spin diffusion. Figure |3ja) 
shows S z (y) for different t with B x = +1T. The ex- 
perimental data (symbols) are fit to the product of 
a Gaussian, Aexp(— y 2 /4w 2 ), and a cosine function, 
cos (27r(2/ — 2/o)Mso) (solid lines). The Gaussian enve- 
lope describes the solution of the spin diffusion equation, 
for which w 2 = D s t (D s is the spin diffusion constant). 
Spin dephasing is included in the time dependence of the 



amplitude A(t). 

A linear fit of w 2 to D s t provides a direct measure of 
D s [Fig. ||;b)], and we obtain D s = (385 ± 15) crmV 1 . 
Note that our experiment uses the spin as a label of the 
electrons and therefore tracks the spin and not the charge 
diffusion. The sensitivity of spin diffusion to electron- 
electron interactions p~4j[15] explains the ten- fold smaller 
value of D s as compared to the charge diffusion constant 
D « 4000 cm 2 s -1 (as calculated for the measured elec- 
tron mobility \i ~ 22m 2 /Vs and a sheet carrier density 
n s = 5 x 10 15 m -2 ). 

Figure (3^b) shows the decay of A(t) with t. The dif- 
fusive expansion of the excited spins into an area pro- 
portional to D s t decreases S z proportional to (Dst) -1 . 
Additional spin decay is mainly induced by deviations 
from perfect SU(2) symmetry [16 j and can be described 
by an exponential decay proportional to exp (— t/r 3 ), 
where r s is the spin lifetime. A(t) is therefore fit to 
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FIG. 3: Spin diffusion and SOI characterization, (a), Normalized line scans S z (y) for various t. The experimental 
data (symbols) are fit to Aexp(— y 2 /Aw 2 ) x cos (2n(y — yo)/\so) (solid lines). The parabola follows ±2\/D s t. (b), Symbols: 
Amplitude A and squared width w 2 of the Gaussian envelope. Solid lines: fits to (t — to) 1 exp(— t/r s ) and D s £, respectively. 
Error bars: Uncertainty of w 2 in a single fit. (c), Shift yo of the spatial oscillation for —1, and 1 T. (d), Symbols: PSH period 
Aso for and 1 T. Solid line: smoothed data at 1 T. (x) and (+): fits to data at three- and nine- fold lower pump intensities. 



(t — to) -1 exp(— £/r s ), yielding r s = (1.1 ±0.1) ns. This is 
about 30 times longer than the Dyakonov-Perel spin de- 
phasing time calculated from the measured D s and SOI 
strength, in agreement with a spin decay time r z « 35 ps 
obtained for 20-/im-wide laser spots, where the spatial 
correlations of the PSH are averaged out (see supple- 
mentary information). 

The measured r s is limited by the two SU(2)-breaking 
contributions, namely, the cubic Dresselhaus SOI Ps and 
the imbalanced SOI \a\ ^ \Pi — Ps\- It can be shown that 
(supplementary information) 

2 

T7 1 « 2D s ^(Wi + (a -/?i+/3 3 ) 2 ), (1) 

where m = 6.1 x 10 -32 kg is the effective electron mass. 
From the measured r s and D S1 we find 3/3f + (a — j3\ + 
/3 3 ) 2 = 1.7 x 10 -26 eV 2 m 2 . This relation restricts /3 3 to an 
upper limit of 0.77 x 10 -13 eVm, which would be reached 
for a = Pi— Ps- From the measured r s alone, it is not pos- 
sible to differentiate cubic Dresselhaus contributions from 



imbalanced SOI. However, as we will discuss later, the 
latter can be separately determined from a small asym- 
metry in S z (x, t) maps that appears if B is applied along 
y. This will allow us to quantitatively describe all SO 
coefficients. 

The fitted PSH period A SO is shown in Fig.gd). A de- 
crease from & 10 jam at t = 200 ps to ~ 7.3 jam at 1.5 ns is 
well represented in several periods of S z (y) and must be 
related to a continuous change of Aso with t. From the 
relation Aso = 7rh 2 m~ 1 (a+f3i— Z^) -1 , we determine that 
\a + /3i — Ps\ increases from 3.5 to 4.9xlO _13 eVm during 
this time. Starting with a broad positive peak in S z (y) 
at t = 10 ps, the spin helix continuously adapts to the 
decreasing Aso- The initially weaker SOI is most likely 
induced by the photo-excited charge carriers that recom- 
bine with time, and is possibly also affected by cooling 
of hot electrons. The SO coefficients are sensitive to the 
screening of the confinement potential and modifications 
of the Fermi energy: Both an increase of Ps with charge 
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density and a reduction of a and ft with screening could 
explain the observed decrease of Aso with t. Supporting 
this interpretation, measurements at lower pump inten- 
sities yield an initially smaller Aso [Fig. [3f d)])- 

Aso is found to be independent of B x [Fig. [3jd)]. To- 
gether with the insensitivity of A(t) on B x , this demon- 
strates the decoupled influence of the Zeeman and the SO 
energy on the electron spins in the case where a magnetic 
field is applied along the unidirectional Bso- 
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FIG. 4: Dependence of the total magnetic field Btot 
on k. (a), The x-component of BgQ is shown (arrows) as 
a function of k for a > 0, Pi — ft > and g < 0. (b), 
If an external magnetic field B is applied along x, \Btot\ is 
different for k y > and k y < 0. The arrows represent the 
size and direction of B^q x + B x . (c), For B || y, \B to t \ is the 
same for ±k y . (d), Away from the SU(2) symmetry point 
(a 7^ ft — ft), the ^/-component of BgQ is non-zero, and the 
superposition with B y leads to a different £?tot for k x > 
than for k x < 0. This allows us to determine a — Pi + ft. 



To understand this decoupling better, it is instruc- 
tive to plot the directional dependence of Bso on k = 
k(cos sin#). For (OOl)-grown QWs, Bso is in the (x,y)- 
plane for all k. It can be written as the sum of two 
terms, BgQ and BgQ. The former is responsible for 
the PSH formation, whereas the latter leads to spin de- 
phasing (see supplementary information). The x- and 
^/-components of BgQ are proportional to k y (a -f ft — ft 3) 
and —k x (a — Pi + ft), respectively. In our case, \a + 
Pi — Pz\ ^> \ol — Pi + ft|, which means that B$q x drives 

the PSH [Fig. [4^ a)], whereas a remaining B$q y breaks 
the SU(2) symmetry and leads to spin dephasing. In the 
following discussion, we omit the superscript and mean 
BgQ when we write Bso- 

As shown in Fig. [4^b), B x < superposes with Bso, x 



such that the total field, B tot — |B + Bso I? is larger for 
k y > than for k y < 0. Translating the momentum hk y 
into a position y using y = hk y t/m, this is exactly what 
is seen in the measurement at B x — —IT in Fig. |2|d): 
The total precession frequency v of an electron with mo- 
mentum hky is determined by the sum of the Zeeman 
splitting, g[i^B xl and the SO splitting, 2k y (a + Pi — P%). 
Electrons with k y = remain at y = and precess with 
v — \g^BB x /h\. There is a path yo(t) [dashed line in 
Fig. [2jd)] on which Zeeman and SO energies cancel each 
other and consequently the spins do not precess. This 
path is characterized by 



dyo/dt = -hgfi B B x /[2m(a + ft - ft)]. 



(2) 



In agreement with this equation, the fitted yo in- 
creases linearly with t and dyo/dt changes sign with B x 
[Fig.j3[c)]. Inserting the measured dyo/dt and g into 
Eq. ([2]), we obtain a + ft - ft = 4.8 x 10" 13 eVm, con- 
sistent with the value obtained from Aso, but here the 
sign is directly determined by the sign of dyo/dt. 

We now investigate how close the SOI in our QW is 
tuned to the SU(2) symmetry point. For that purpose 
we apply B along y. Figure [4jc) shows that then, B tot 
is not unidirectional anymore. Even though a map of 
S z (y,t) at By = -IT [Fig. [5ja)] still indicates that a 
helical mode evolves for t < 400 ps, S z quickly decays 
at longer t. In addition, an oscillation with t is seen at 
(x,y) = (0,0). This is related to the precession of S 
about By. On the other hand, the superposition of B y 
and Bso, y [Fig.Qd)] leads to an asymmetry for opposite 
signs of k xi which can be observed in maps of S z (x,t) 
[Fig. |5jb)]. We find a small tilt of the positions xo(t) of 
constant spin precession phase. In analogy to Eq. 
the tilt is given by 



dxo/dt = hgnBB y /[2m(a - ft + ft)]. 



(3) 



dxo/dt is a measure of the detuning from balanced SOI. 
We obtain dxo/dt w — 280/im/ps for B y = —IT, and 
about twice that value for B y = — 2 T [dashed lines in 
Fig. [5^b)] . From this, it follows that a — ft + ft « 
—0.3 x 10 _13 eVm, where the sign is directly determined 
by the sign of dxo/dt. 

We can now derive the size of all SO coefficients in our 
sample. From the known sum and difference of a and 
ft — ft, respectively, we find a = (1.6 — 2.3) x 10 -13 eVm 
and ft -ft = (1.9-2.6) x lO" 13 eVm. The cubic Dressel- 
haus SOI ft is estimated from Eq. ([T]) based on the ob- 
served r s , yielding ft w 0.7 x 10 -13 eVm. We then obtain 
ft w (2.6 — 3.3) x 10 -13 eVm, which is in excellent agree- 
ment with results from a similar QW structure [6] . Using 
the theoretical expression ft = —j(k^) [J, with (k%) = 
3.7xl0 16 m -2 as obtained by solving the one-dimensional 
Poisson and Schrodinger equations, we determine the 
Dresselhaus coupling parameter 7 w — 9eVA 3 . This is 
in agreement with previous work where 7 was found 
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FIG. 5: Interplay of the PSH with an external magnetic field, (a), Experimentally measured and numerically 
simulated maps of S z (y,t) for B y = —IT. The formation of a helical spin mode is challenged by simultaneous spin precession 
about By. (b), Maps of S z (x,t) recorded at B y = — 1 and — 2T show the spin precession about B y as well as an asymmetry 
of the precession phase with ±x (dashed lines), which is attributed to a slight detuning from the SU(2) symmetry point. 



to be in the range of —4 to — 8eVA 3 [5]. Inserting 
7 -9eVA 3 into f3 3 -\^irn s , with n s 5x 10 15 m~ 2 , 
we find /3 s = 0.7 x 10 _13 eVm, exactly the same value 
as obtained from r 3 . This indicates that the PSH de- 
cay is well described by Eq. ([TJ. The Rashba coeffi- 
cient a can be related to the electric field £?qw m the 
QW by a = tqw^qw, defining a proportionality con- 
stant a = tqw From the calculated conduction-band 
profile, we estimate £?qw = 4 — 6 x 10 6 V/m. Using 
a w 2.3 x 10 _13 eVm, we obtain tqw = 4 — 6eA 2 , in 
good agreement with theoretical prediction pQ and ex- 
periment [6]. 

To crosscheck the consistency of our experimentally 
determined SO coefficients, we simulate spin diffusion in 
the presence of SOI using a Monte Carlo approach that 
combines semiclassical spin dynamics and diffusion (see 
numerical methods). The simulation uses the experimen- 
tally determined values for a±, /3i, j3s and D S1 and repro- 
duces the spin dynamics in the combined field of B and 
Bso remarkably well [see Figs. [5|a) and (b)]. It also 
confirms the analytical relation of Eq. ([2| [dashed lines 
inFig.gb)]. 

Finally we want to discuss the transition from a PSH- 
dominated regime to a regime in which the external 
magnetic field dominates. Figure [6ja) shows two maps 
of S z (y,t) recorded at B y = -0.375 and -2T. At 



By = — 0.375 T, the formation of a helical mode is read- 
ily identifiable, but r s is significantly lower than for 
0T [see Fig. [2fc)]. At -2T, no signature of a heli- 
cal mode is observed and spins precess about B y with 
v = \g/jiBBy/h\. To determine v and r s versus B yi the 
time traces S z (t) at (x,y) = (0,0) [Fig. [6]^b)] are fit to 
(t— to) -1 exp(— t/r s ) cos(27rz4). Starting from the PSH at 
0T, t s rapidly decreases with increasing \B y \ [Fig. [6^c)]. 
Interestingly, spin precession about B is suppressed in 
the regime of \B y \ < 0.5 T. Even though at \B y \ = 0.5 T 
the Zeeman spin splitting (5 fieV) is more than one order 
of magnitude smaller than the SO spin splitting (up to 
170 fieV at the Fermi energy), the symmetry-breaking 
induced by B y dramatically decreases r s from 1 ns to 
below 100 ps [Fig. [6jc)] . In an intermediate regime of 
0.5 T< < 1.5 T, the helical spin mode coexists with 
spin precession about B, and r s recovers from its mini- 
mum value. For > 1.5 T, r s monotonically decreases. 
In the regime where the Zeeman energy approaches the 
SOI energy, we expect that r s approaches the value given 
by the Dyakonov-Perel expression for homogenous po- 
larization. Also shown in Fig. [6jc) are fits to numerical 
simulations of S z (t) at (x,y) = (0,0) that exhibit the 
same features as the measured data. 
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FIG. 6: Detuning from the PSH regime, (a), Left: 
Maps of S z (y,t) for B y = -0.375 T. The PSH is still visible, 
but the spin lifetime has decreased to a few 100 ps. Right: 
Maps of S z (y,t) for B y — — 2T. Instead of the formation of 
a PSH, a spin precession about B y is seen, (b), Time traces 
S z (t) exhibit no spin precession for \B y \ < 0.375 mT, and r s 
decreases rapidly with \B y \. (c), Precession frequency v and 
spin lifetime r s versus B y . Open symbols: Experimental data. 
Filled symbols: Numerical simulation. 



duction band of the QW. The spin polarization compo- 
nent S z along the z || [001] direction is measured by a 
probe pulse at A = 799 — 801 nm using the polar magneto- 
optical Kerr effect [T7HT9] . Both pulses are directed 
through a single high-numeric-aperture lens placed inside 
the cryostat. The spatial position of the pump beam on 
the sample surface is scanned by controlling the angle of 
the incident pump beam with a Galvo mirror. In spatial 
maps of the spin distribution S z {x 1 y) 1 the coordinates x 
and y indicate the position of the probe beam relative to 
the scanning pump beam. A non-linearity in the angle 
control was corrected to obtain calibrated values for x 
and y. The pump pulse is spectrally removed with a low- 
pass filter prior to detection. The focal distance is opti- 
mized to achieve a maximal spatial resolution of 2 /mi in 
diameter. The probe pulse is synchronized to the pump 
pulse and delayed with a mechanical stage by a time t to 
monitor the time evolution of the spin polarization S z (t) . 
The pulse lengths of the pump and the probe beam are 
« 60 ps and 3ps, respectively. Typical power intensi- 
ties of pump and probe were 250 and 50 /iW, respec- 
tively, with pulses arriving at a repetition rate of 80 MHz. 
Fig. |3jd) includes data with the same probe power, but 
100 and 30 /iW pump power. The data of Fig.[6]have been 
taken at a pump and probe power of 150 and 50 /iW, re- 
spectively. All measurements were carried out at a sam- 
ple temperature of 40 K. Photo-excited electrons equili- 
brate rapidly (< lps) to a Fermi distribution, however 
the initial electron temperature may well exceed the lat- 
tice temperature. On a time-scale smaller than 100 ps, 
it is expected that these electrons are cooled by phonon 
emission [20] and by scattering with the cold background 
electron gas. Transport measurements on etched mesa 
structures with Ohmic contacts yielded n s = 5 x 10 15 m -2 
and fi = 22m 2 /Vs. From n s we estimate a Fermi energy 
of 18meV above the QW ground state. The QW ab- 
sorbtion edge lies at 1.54 eV at 40K. In the experiment, 
we have to account for a finite laser spot size, leading to 
a finite width of the spin distribution at t = 0. There- 
fore we introduced a time offset to in the fit function 
(t — to) -1 exp(— t/r s ). This adds an uncertainty in the r s 
extracted [error bars in Fig. [6Fc)]. 



Methods 
Experimental methods 

The 2DEG investigated is confined to a 12-nm-wide 
(OOl)-oriented GaAs/AlGaAs QW placed 95 nm below 
the surface. Asymmetric Si modulation doping provides 
a sufficiently strong Rashba SOI to nearly balance the 
Dresselhaus SOI contribution. We employ a pump-probe 
approach in which a circularly polarized pulse of a mode- 
locked laser at A = 785 nm (full- width at half-maximum 
of 10 nm) excites spin-polarized electrons into the con- 



Numerical methods 

A numerical simulation of spin dynamics is performed 
using a two-dimensional Monte Carlo approach in which 
at time t = 0, 20000 spins are distributed at coordi- 
nates (x, y) with a Gaussian probability distribution of 
width 2 /mi. Wave numbers (k x ,k y ) are uniformly dis- 
tributed on the Fermi disc, and the spins S are oriented 
along the z-direction. In discrete time steps, updates 
of (x, y) and S are calculated, treating the spin preces- 
sion about the sum of Bso + B semiclassically. Scat- 
tering is accounted for by isotropically redistributing the 



8 



charge carriers on the Fermi disc with scattering prob- 
ability r = 2D s /vp, where is the Fermi velocity. 
For the simulation, we used a single set of SO coeffi- 
cients, a = 1.7, Pi = 2.7, ^ 3 = 0.7 x 10" 13 eVm, and 
D s = 385 cm 2 s _1 to reproduce the measured S z as shown 
in Figs. [5j [6] and [8j The variation of the SO coefficients 
with t, as observed in the measurement, has not been 
included in the simulations. 



The relaxation time r s of the helical mode can be ex- 
pressed [16] in terms of squared components (B^ Qx ), 
(^io,y)> wnere (•••) denotes averaging across the Fermi 
disc, as well as (Bso yX k y ). Close to the SU(2) symmetry 
point, i.e., for a w Pi — Ps, the contribution from BgQ to 
t,T 1 is proportional to (a — Pi + ^3) 2 and the one from 
to 3/3| [see Eq. (1)]. 



Supplementary information 

The effective spin-orbit magnetic field and 
dephasing of the persistent spin helix. 

The spin-orbit interaction for a two-dimensional elec- 
tron gas confined in a (OOl)-grown GaAs QW is de- 
scribed by the three spin-orbit parameters a, Pi, and 
Ps- It can be expressed as an effective spin-orbit mag- 
netic field Bso, which in the coordinate system x || [110] 
and 2/ || [110] is given by [21] 




300 



t(ps) 



>so 




with Pi = —j(k 2 } and Ps — —jk 2 /4. {k 2 ) is the expec- 
tation value of k 2 with respect to the QW ground-state 
envelope wave- function, k = \j2i\n s is the Fermi wave 
number of the 2DEG with sheet density n s , and 7 is the 
Dresselhaus coupling parameter. It is convenient to write 



FIG. 7: Dyakonov-Perel dephasing rate Time trace 
( 4 ) S z (t) of spin polarization for a measurement with 20-/im- 
wide laser spots. In this case, spatial spin-orbit correlations 
are averaged out and S z (t) decays with the Dyakonov-Perel 
dephasing rate. 



> SO 



>(3) 
^SO 



E21, with 



B (i) _ 
^so - 



2k 
9Vb 



(a + P ± - P 3 )sm0 
-(a - P\ + p 3 ) cos( 



and 



Dyakonov-Perel dephasing rate 



In the diffusive limit, where the scattering length is 
much smaller than the spin-orbit length and if spins per- 
(5) form a random walk on the Bloch sphere (Dyakonov Perel 
regime), the expected decay rate for a spin polarization 
along the z-direction is [21] 



R (3)__2W fosinSe \ 

Here is the angle between k and the x-axis. Cross- 
correlations (Bso^ x k y ) and (Bso^ y k x ) between compo- 
nents of Bso and components of k are responsible for 
the formation of helical spin modes. As BgQ disappears 
in these terms, the spatial pattern of the spin modes can 
be described by BgQ alone. In general, a PSH forms if 
a = ±(/?i — Ps). In our sample, BgQ mainly points along 
the ^-direction, which according to Eq. ([5| is the case if 
a and Pi — Ps have the same sign. In this situation, the 
helicity of the mode is determined by the sign of the cor- 
relation (B S o,x k y ) = k 2 /(gii B ) ■ (a + Pi - Pz)- From the 
shift of the PSH in an external magnetic field [Fig. 2(d)], 
we determine a positive helicity, i.e., a uo + mode, and 
therefore a + pi — Ps > 0. 



T" 1 = 8D s m 2 h- 4 [a 2 + (ft - fof + ft] . (7) 

If the spatial correlations of the PSH are averaged out, 
r s is given by Eq. ([7]) instead of Eq. (1). Using 20- /im- 
wide laser spots, we measure r z « 35 ps at 30 K (Fig. [7]). 
This is in perfect agreement with the value obtained from 
Eq. using a = 2.3 x 10- 13 eVm, Pi - p 3 = 2.6 x 
10- lJ eVm, p 3 = 0.7x 10- 13 eVm and D s = 385 cmV 1 . 



Numerical simulation data 

To compare experiment and simulation, we present the 
numerically simulated S z (y, £)-maps for B || y that cor- 
respond to the experimental data shown in Fig. [6] Simu- 
lation and experiment agree remarkably well in both the 
Bso- and the B-dominated regime (Fig. [§]). 
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FIG. 8: Detuning from the PSH regime. Map of S z (y, t) 

for By — —0.375 and — 2T. Top panels: Experimental data 
(same data as shown in Fig. 6). Bottom panels: Numerical 
simulation. 
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